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Abstract. The non-linear evolution of a magnetic island is studied using the 
Vlasov gyro-kinetic code GKW. The interaction of electromagnetic turbulence 
with a self-consistently growing magnetic island, generated by a tearing unstable 
> 0 current profile, is considered. The turbulence is able to seed the magnetic 
island and bypass the linear growth phase by generating structures that are 
approximately an ion gyro-radius in width. The non-linear evolution of the 
island width and its rotation frequency, after this seeding phase, is found to 
be modified and is dependent on the value of the plasma beta and equilibrium 
pressure gradients. At low values of beta the island evolves largely independent 
of the turbulence, while at higher values the interaction has a dramatic effect on 
island growth, causing the island to grow exponentially at the growth rate of its 
linear phase, even though the island is larger than linear theory validity. The 
turbulence forces the island to rotate in the ion-diamagnetic direction as opposed 
to the electron diamagnetic direction in which it rotates when no turbulence is 
present. In addition, it is found that the mode rotation slows as the island grows 
in size. 


1. Introduction. 

Magnetic islands in a tokamak can lead to loss of confinement or even major 
disruptions of plasma. The tearing mode mm, specifically the neoclassical tearing 
mode (NTM) [3] is expected to set the beta limit in a reactor [H 0 IS] ■ 

The tearing mode is centred around flux surfaces with rational values of the safety 
factor. Current perturbations at these surfaces, coupled with non-ideal effects such 
as electron inertia or resistivity, can become unstable. This process, via magnetic 
reconnection, forms a growing magnetic island with its corresponding radial magnetic 
field component. When collisions can be neglected, the mechanism of reconnection is 
through the electron inertia, where the singular layer width (the layer where magnetic 
reconnection can occur), and in turn the growth rate, are related closely to the electron 
skin depth [7], this corresponds to the collision-less tearing mode [8]. When collisions 
become significant, it is plasma resistivity which determines the singular layer width. 
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producing a classical collisional tearing [T] or, more relevant for present day tokamaks, 
a semi-collisional mode [Hiiiiiin]. In high temperature, weakly collisional, plasmas the 
mode grows on a slow time-scale. 

The stability and evolution of magnetic islands is the result of the interplay of 
variety of processes and their interaction with electromagnetic turbulence is further 
expected to modify their stability and evolution Turbulence also generates, and 
is regulated by zonal and meso-scale flows [12] , a process which is further complicated 
by the presence of magnetic islands [m El da dg Hz]. Plasma turbulence, zonal 
flows and tearing modes occupy disparate time and length scales, however, early in 
their evolution magnetic islands can be very narrow and thus comparable to turbulent 
length scales, as such, their evolution can not be considered to be independent of 
the turbulence EldH], and has been shown to be influenced by turbulence in both 
gyro-fluid simulations [iniHO] and gyro-kinetic simulations in toroidal geometry pTj . 

Turbulence has been suggested as a source of seed islands for neoclassical tearing 
modes which display no clear threshold behaviour (i.e they only grow after an initial 
perturbation is formed) [22], and can also be triggered by an island generated by an 
unstable classical tearing mode [23] . 

The large separation in time-scales, overlap of length scales and the complexity 
of the multitude of interactions makes analytical theory difficult, therefore here we 
approach the problem using massively parallel, state-of-the-art kinetic simulation. 
The rest of this paper is organised as follows. In Sec. [2] we outline the simulation set¬ 
up and the parameters considered. In Sec. [3] the non-linear evolution of the tearing 
mode is described without turbulence present. In Sec. |4]the non-linear growth of the 
island self-consistently interacting with electromagnetic turbulence at various values of 
/3e is considered (/3e is the ratio of the plasma pressure to the magnetic pressure). The 
effects of turbulence on the island rotation frequency is considered in Sec. [5] Finally, 
the results are summarised in Sec. n 


2. Parameters and set-up. 


The global version of the electromagnetic non-linear gyro-kinetic code GKW [24] was 
utilised in this study. A full description of the equations solved, and the implementation of 
the tearing mode current drive can be found in references [251 HI]. 

For all the simulations presented here the normalised gyro-radius p* = pi/R — 0.005 
is used. This gives a n = 1 toroidal mode number with normalised binormal wave-number, 
k(^pi = 0.053. In a full tearing mode and turbulence simulation a total of 36 toroidal modes 
were used giving a maximum kc^pi = 1.86. The ion and electron temperatures are assumed 
to be equal, Ti = Te. Resolutions in the parallel, parallel velocity, magnetic moment, 
radial directions are, Ns = 64, = 64, Nfj, = 16, Nx = 256 respectively. This radial 

resolution places at least three radial grid points per singular layer width depending on 
the value of /3e. The linearly growing tearing mode is found to be a semi-collisional mode 
[8]. As such, the singular layer width is closely related to the electron skin depth, which 
is given by, — = —^ ^ ^ which must be resolved. The /3e is the electron /3 defined as 


Pe = neTsj{BQl2po), where and are equilibrium density and temperature and Bq is 
the magnetic field strength at the outboard mid-plane magnetic axis. 

The tearing mode instability is driven by a non-homogeneous current density profile. In 
our implementation this is introduced by applying an electron flow profile in the equilibrium 
which is calculated to be self-consistent with the q-profile. The current proHle used in this 
paper is the same model as used by Wesson et. al. | 26lI27| . Here the current density profile 
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Figure 1. The radial profiles of the (top) linear Aj| and (j) eigenfunctions at 
/3e = 0.05% and the A|| profile at non-linear saturation (dashed line). (Upper 
middle) the safety factor (q) and magnetic shear (s) profiles. (Bottom middle) 
the (grey solid) equilibrium electron density profile (black) the time averaged 
electron density profile through the island 0-point and (dashed) a snapshot 
electron density profile. (Bottom) the perturbed electron current profile at island 
saturation. (Blue) vertical dashed lines represent the approximate location of the 
island seperatrix. 


has a parabolic form, j = 20(1 -\- {r/aY)^ which introduces an electron flow which enters the 
equations via the source term, 

NFmg = —2—^VUe-FMe, (l) 

where Me = Me(r) is the electron flow which is related to the current density profile by, 
j = n^u^e. In the rest of this paper, when a simulation without the drive is mentioned, it is 
terms of this form of Eq. [T] that are neglected from the equations solved. The safety factor 
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profile is calculated from the current profile, which has the analytic form, 
Bs r'^/c? 


= (la 


R 1 — (1 — r'^/a'^Y 


( 2 ) 


where qa is the value of the safety factor at radial co-ordinate ip = a/R. An example of the 
safety-factor and magnetic shear (s = rjqdqjdr) used can be seen in the middle panel of 
Fig. [T] In the simulations presented here, q, at the at the plasma edge is set to qa = 3.5 
and the safety on the axis is 5 ( 0 ) = 1.75. Therefore, q = 2, q = 2.5, q ~ 3 and many more 
higher mode number rational surfaces {q = m/n, where m is the poloidal mode number) 
exist within our computational domain. The radial domain covers 0.1a < r < a where a is 
the minor radius. 

The safety factor profile given above is found to provide a linearly unstable semi- 
collisional tearing mode |25| . with a positive stability parameter [T] defined as. 


A' = 


1 dA« 


An dr 


(3) 


across the singular layer, whose position is denoted by rs and upper and lower boundaries by 
and . This parameter represents the discontinuity in dA\\/dr and measures whether 
mode growth is energetically favourable (A' > 0 for the mode to grow). The tearing mode 
has the familiar [2H] linear structure of which is shown in Fig. [T] showing the position of the 
resistive layers at the q = 2 and q = 3 rational surfaces. 

While a full Fokker-Planck collision operator is implemented in GKW, here we use only 
pitch-angle scattering of the electrons off the bulk ions for simplicity. The normalised collision 
frequency (normalised to the trapping/de-trapping rate), = 4z^ei /3yJ~{ne^) = 0.12, where 

the collision frequency is defined as a* ■ The density and temperature profiles 

47r€Q mg Ug 

have the radial form. 
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where R/Ln = —{R/n)dn/dr and R/Lt = —{R/T)dT/dr are the logarithmic gradients, 
Ar is the length scale on which the gradient increases from zero to this at the boundaries. 
This form maintains a constant gradient across most of the radial domain, while dropping 
the gradients to zero close to the boundaries to prevent numerical issues. In all cases the 
simulation was initialised with grid-scale noise at a small amplitude. 


3. Non-linear evolution and saturation 


Figure [2] shows (black dashed) the time evolution of the magnetic island width for a non¬ 
linear simulation where the background pressure gradient is insufficiently large to drive the 
ITG mode unstable. The magnetic island width is related to the magnitude of the vector 
potential (in the constant-?/) approximation) by the expression 


w 



(4) 


where Ayo is the mode amplitude at the singular surface. 

The parameters utilised were R/Ln = 1.0, R/Lt = 3.5 and /3e = 0.1%. In this case 
the island evolves in the manner as predicted by Rutherford [^. At small time, the island 
grows exponentially at its linear growth rate until it reaches a width (at ~ SOOR/vthi) that is 
approximately the singular layer width (~ 0.7pi Parameters are shown in Table[T] Case 2.). 
After this the growth of the island slows down and enters the non-linear regime |29| where the 
growth of the island becomes algebraic rather then exponential {w oc t). Eventually the island 
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Figure 2. (left) The time evolution of parallel vector potential amplitude (A||) 
at the q = 2 resonant surface for (black dashed) /3e = 0.1% with R/Ln = 1.0 and 
R/Lj- = 3.5, (blue) /Jg = 0.05% with R/L„ = 1.5 and RjLx = 5.0 the dashed 
line represents the linear scaling of the unstable TM with these parameters. The 
solid black line is the same simulation as the dashed black line, but transposed 
so that the amplitudes match at the end of the island seeding phase, (right) A 
snapshot of (top) (p and (bottom) A|| at island saturation. Vertical dashed lines 
correspond to the radial positions of the singular layers. 


reaches a saturated width which is calculated here to be 3.7pi, consistent with quasilinear 
calculations for this form of equilibrium current profile [30]. Also shown in Fig. [T] is the 
electron density profile, where the flattening of the profile through the O-point is evident at 
both the q = 2 (rnjn = 2/1 mode) surface and the smaller island at the g = 3 (mjn = 3/1) 
surface. 

Stochastisation of magnetic islands is a well known phenomenon, and it has also been 
shown in non-linear gyro-kinetic simulations of magnetic islands with turbulence that higher 
mode numbers may excite and cause stochastisation of the island seperatrix and an almost 
complete break down of the island X-points [21]. Here however it can be seen that even when 
the equilibrium gradients are too small for turbulence to form, the non-linear modification 
of the density and temperature profiles can cause a steepening of gradients near the island 
seperatrix (also shown in Fig. [T]) and this, in turn, can cause a high-k ITG to become 
unstable in the island vicinity. Figure [21 shows four time slices of the island electro-static 
potential during the non-linear phase of the island evolution. It is evident that a high mode 
number, drift-wave becomes unstable at the island seperatrix. A similar process has been 
observed in 2D fluid simulations of an tearing mode m- These small scale ITG modes can 
then interact and form turbulence, which in turn can cause the island seperatrix to become 
stochastised. The bottom right panel shows a Poincare plot of the magnetic field at the 
saturated state, clearly showing a stochastic layer around the island which, initially, had a 
coherent form (bottom left). Also visible are islands generated at higher order mode numbers 
(e.g. m/n=8/3, 9/4) that are linearly stable, but here are non-linearly excited and causing 
stochastisation from island overlap |32[I33 |. 
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Case 

fie (%) 

R/Ln 

RILt 

Growth rate, 

7 (vth/R) 

Frequency, 

UJ (Vth/R) 

Electron 
skin depth 

iP^) 

1 

0.05 

1.5 

5.0 

0.045 

-0.16 

1.04 

2 

0.1 

1.0 

3.5 

0.032 

-0.12 

0.73 

3 

0.1 

1.5 

5.0 

0.024 

-0.15 

0.73 

4 

0.2 

1.5 

5.0 

0.0090 

-0.15 

0.52 

5 

0.1 

2.2 

6.9 

0.0032 

-0.28 

0.73 


Table 1. Table of the growth rates, mode frequencies and the corresponding 
electron skin-depths for the corresponding parameters in simulations utilised in 
this paper. All linear calculations are for the n = 1 toroidal mode, and the most 
unstable mode is always found to be a tearing mode. 


4. Turbulence interaction 

As the equilibrium density and temperature gradients are increased, the ITG becomes 
unstable and electromagnetic turbulence forms. A snapshot of a. 41 — (^ plane during a 
turbulence simulation is shown in the right hand panels of Fig. [21 They show the electro-static 
potential (top) and electromagnetic potential (bottom) at point where the island saturates 
t ~ 700R/vthi. The turbulence, in turn, interacts with a growing magnetic island. 

Firstly we concentrate on a simulation with an electron j3e = 0.05%, (Case 1 in Table[T]) 
here the turbulence is fundamentally electro-static, and the tearing mode has a large singular 
layer width (Since he ~ ~ pi). The equilibrium temperature and density gradients 

have the logarithmic gradients, R/Lt ~ 5.0 and R/Ln = 1.5 respectively. Gradients are at 
the reference radius, r = ro = 0.5a and are the same for both ions and electrons. These 
gradients are large enough to drive micro-instabilities which lead to turbulence but not large 
enough to completely stabilise the linear tearing mode [341 125] . an observation that was 
confirmed by linear simulations of the n = 1 toroidal mode using the same parameters. The 
calculated growth rates and frequencies are shown in Tabled] 

The time trace of the n = 1 toroidal component of the electro-magnetic potential (A||) 
for these simulations are shown in figure [2| Also included is the trace of the n = 1 electro¬ 
static potential, f (Green line) showing that it is largely determined by the turbulence, 
reaching rapidly a saturated state. In contrast, the island is still growing as evidenced by 
the increasing A||, the electrostatic potential, has saturated. It can be seen comparing 
the trace with turbulence (blue line) to that without (black dashed line), that the turbulence 
provides an initial island structure |20l 1351136 | that is approximately the size of the ion-gyro- 
radius. A gyro-kinetic study of this seeding mechanism can be found in |21| . Without a 
background current profile (grey line) the amplitude of the mode saturates and maintains 
this amplitude, but with a tearing unstable current profile the seed island structure starts to 
grow (from approximately t = lOOR/vthi) and eventually saturates at a half width that is 
calculated to be ic = O.Tpi after JOOR/vthi- 

Taking the curve of the island evolution without turbulence present, and shifting the 
x-axis so that the mode amplitudes are comparable at the end of the seeding phase (approx. 
lOOR/vthi). It can be seen that the non-linear island evolution after this point agree almost 
exactly from this point until island saturation. The saturation amplitudes also have a 
very close agreement. While the equilibrium pressure gradients in these two simulations 
are slightly different, giving different linear growth rates, this is not expected to change 
significantly the non-linear behaviour of the island, therefore a comparison can be made. It 
is seen that the turbulence seeding mechanism bypasses the linear regime entirely, providing 
a structure much bigger than the singular layer width in almost all the cases considered, and 
approximately the layer width in this low beta case. Above this width the island evolves. 
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Figure 3. Snapshots of the non-linear evolution of the electro-static potential of 
a growing magnetic island when the equilibrium gradients are insufficient to drive 
micro-instabilities. However, the growing island can induce a small scale drift- 
wave around its separatrix, causing a stochastisation around the seperatrix. This 
stochastisation can be seen in the Poincare plot (bottom right) which corresponds 
to the last time slice at island saturation, (bottom left) is the same but during 
the linear phase of the magnetic island evolution. The q = 2, q = 2.5 and q = 3 
resonant surfaces are also plotted and higher order islands are also visible. 


with this value of /3e, almost independent of the turbulence around it and is governed by the 
positive stability parameter A'. 

However, as Pe is increased, a marked change in the island evolution is observed. Figure 
[4] shows the time evolution of the island width in simulations where Pe was increased to 
0.1 — 0.2% (Cases 3, 4 and 5 in Table]!}. At early time a similar evolution to the low Pf. case 
is seen. The electromagnetic component of the turbulence provides a seed magnetic island 
structure which further grows, once turbulence has saturated, if there is a tearing unstable 
current profile present. However, after this point a different evolution is seen. Without 
turbulence the island evolves non-linearly with an algebraic growth, w oc t. With turbulence, 
in all three simulations presented here varying Ps and equilibrium pressure gradients, the 
island grows exponentially at a rate that is the same as its corresponding linear growth rate. 

The linear growth scaling of these simulations are shown as corresponding colour dashed 
lines in the hgure. The calculated growth rates and frequencies are outlined in Table [T] 
Two stabilisation effects are evident. Firstly, as the density and temperature gradients are 
increased (for constant /3e) the linear growth rate drops due to the generation of diamagnetic 
currents outside the singular layer, which cause a stabilisation of the tearing mode |34| . 
Pressure gradients, combined with magnetic curvature are also known to have a stabilising 
effect on the linear m and non-linear mode growth [SS], however at these low P^, the effects 
are likely to be small. Secondly, as the Ps is increased, the growth rate reduces due to the 
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Figure 4. The time evolution of the n = 1 parallel vector potential A|| for 
four different simulations where the equilibrium temperature gradient, density 
gradients and electron /3e are varied. Values are (blue) = 0.1%, RILx = 5.0, 
RjLn = 1.5 and (green) /3e = 0.2%, RILj- = 5.0, RjLn = 1-5 and for (magenta) 
/3e = 0.1%, RILt = 6.9, R/Ln = 2.2. (black) uses the same parameters as the 
blue traces but without a current profile. The dashed lines are the corresponding 
linear island growth scalings. Horizonal (black) dashed and (red) dashed lines are 
the island saturated width and singular layer width respectively from a simulation 
with only two toroidal modes (shown in inlay). 


narrowing of the singular layer width. Drake and Lee showed that for a collisionless tearing 
mode (semi-collisional is a slight modification of this) that the growth rate 7 ~ where he 
is the skin depth. A higher /3e reduces the skin depth and thus the growth rate, at lower /3e 
the skin depth is wider and thus a larger growth rate is seen. 

In all cases the electromagnetic turbulence supplies a seed island structure of 
approximately an ion gyro-radius in width, greater than or close to the size of the singular 
layer width (Values are also outlined in Table [!} and thus linear growth scaling should 
not be valid here. The island, when turbulence is not present or in the lower /3e case 
seen in the previous section, at this size grows algebraically. However, the electromagnetic 
component of the turbulence here is larger and thus a stronger interaction between the island, 
turbulence and its corresponding zonal flow is evident [m ESI- This interaction has the 
effect of disrupting or modifying the threshold of the non-linear mechanism which, without 
turbulence, would slow the islands evolution. Thus, due to the interaction of the island with 
turbulence, the growth rate is significantly faster than expected. 

5. Island rotation. 

The amplitude and the sign of the island’s rotation frequency, uj, is of direct consequence 
to models of its non-linear stability |39l 1401141| . The polarisation current [42], one of the 
mechanisms invoked to explain NTM threshold behaviour [43] , is closely related to the mode 
frequency. Whether or not the polarisation current is stabilising is dependent on the sign of 
the frequency | 441145 | . As such, a knowledge of the islands rotation frequency, particularly 
when the width is small, is critical in understanding island stability. Moreover, in the case of 
finite overlap of the ion orbits with the island, the rotation of the island is known to have an 
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Figure 5. (top) Mode amplitude and (bottom) mode frequency for a growing 
magnetic island in the absence of turbulence. The black dashed lines denote the 
approximate amplitude and time where the island starts to evolve non-linearly. 
The red horizontal dashed line denotes the linearly calculated mode frequency. 


effect on the density profile through the magnetic island, which can enhance or counteract 
the flattening of the profile depending on its sign [411146j . This effect has been observed also 
in the presence of turbulence [33 ■ The bootstrap current profile, the drive of the Neoclassical 
tearing mode nisi, is directly related to the density profile. 

When turbulence is not present the magnetic island, in the low-collisionality linear 
regime, rotates at the electron diamagnetic direction [34] and then the rotation frequency 
reduces in amplitude as the island grows and enters the non-linear regime. Since toroidal 
angular momentum is conserved in the absence of external sources, as the island grows the 
plasma mass within the island increases and thus its frequency must decrease [49] ■ This point 
in the simulation is denoted by the grey dashed lines. This is shown in Fig[5]for a simulation 
of an island without turbulence (Case 3 in Table [T]) which shows the mode amplitude (and 
therefore island width) and its frequency as a function of time. 

However, the above description is significantly altered when turbulence is present. Figure 
[^shows the time trace of phase, of the O-point as a function of time for case 1. Where the 
phase is related to the mode frequency by tJ = —dC^/dt. This is achieved by knowing that the 
mode with the longest binormal wave-vector is the most tearing unstable and thus produces 
the dominant island. The position of the O-point is found by Fourier filtering away all higher 
k modes and then finding the minima of the mode in the binormal direction. This location, 
and thus the island phase, can then be traced through a turbulent simulation. For reference, 
an increasing phase with time is in the electron diamagnetic direction, while the opposite is 
in the ion direction. 

Initially, the n = 1 mode has a frequency in the electron direction, as expected for a 
linearly growing tearing mode (rotation at the linear frequency is shown by the red dashed 
line), or micro-tearing mode that can be found at this mode number when the current drive 
term is neglected. Once the ITG grows and turbulence is established it can be seen that 
the island mode begins to rotate in the ion diamagnetic direction. The the black dashed 
line represents the rotation at the ion-diamagnetic frequency (cui* = n = 1 

toroidal mode, and it is evident that the island rotates at approximately this frequency. 
However, it is also visible that when there is an island growing, the mode frequency is also 
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Figure 6. (top) The time evolution of the phase of the magnetic island O- 
point for simulations where electromagnetic turbulence is present in the low fie 
case. The (red dashed) line represents the scaling with the linear frequency. 
(Black dashed) The scaling of the phase with the ion-diamagnetic frequency (oJif). 
(bottom) for reference the amplitude of the n = 1 electromagnetic mode. 


a function of the island’s width. The island slowing down as it grows and reaching almost 
a zero frequency state when it saturates. It may also, as is visible here for the Pe = 0.05% 
case, temporarily change sign but maintaining a small frequency amplitude. 

Figure [7] also shows the time traces of the island phase for cases 3, 4 and for the same 
parameters as case 3 but without the current gradient drive term to form a growing island. 
Here it can be clearly seen that the rotation frequency of the mode when there is no growing 
island is approximately constant, while when a growing island is present in the non-linear 
phase, the frequency is a function of the island width. The frequency decreases in amplitude 
as the island grows. This island deceleration does not seem to change significantly with 
increasing fie or between the different equilibrium gradients between these two simulations. 

6. Conclusions 

Using state-of-the-art gyro-kinetic simulations this paper has presented the non-linear 
evolution of magnetic islands driven unstable by the classical tearing mode, and their self- 
consistent interaction with electromagnetic turbulence. 

A number of simulations were performed varying the plasma /3e and the equilibrium 
temperature and density profiles, but in all cases it was found that the turbulence provides, 
through non-linear interactions, a seed island structure. However, how the island evolves 
after this phase is dependent on the plasma electron fie- At lower beta, the island evolves 
non-linearly largely free of interaction with the turbulence, confirmed by comparison with 
simulations of an island an the absence of turbulence. However, at higher values of beta, the 
interaction becomes stronger and has the effect of causing the island to grow exponentially, 
even though it is no longer in a regime where this is expected because the island widths are 
significantly larger than the singular layer width. 

A magnetic island will rotate in the electron diamagnetic direction when no turbulence 
is present at low collisionality. This rotation slows down as the island grows, through 
conservation of momentum, since a growing island captures a larger plasma mass as it grows. 
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Figure 7. The time evolution of the position of the magnetic island O-point for 
simulations where electromagnetic turbulence is present. The (blue dashed) line 
represents the rotation at the ion diamagnetic velocity for the n = 1 mode. The 
blue trace is for a simulation where the island is not growing as no equilibrium 
current profile is imposed, while the red and black curves represent simulations 
where magnetic islands are growing from a seed, Cases 4 and 3 respectively. 


However it is observed that, when present, the turbulence forces the magnetic island to rotate 
in the ion-diamagnetic direction even when the island is small. When a growing island is 
present, as it grows larger, the rotation slows down, eventually stopping and periodically 
changing sign as it approaches a saturated width. 

Non-linear simulations with turbulence may excite modes that are linearly tearing stable 
and form higher mode number islands. The presence of turbulence and higher order modes 
can cause the separatrix and X-point to become stochastised. However, it was also observed 
that even when the ITG mode is stable, the island’s influence on the background equilibria 
(flattening the proHles within, and steeping just outside) can cause a high mode number ITG 
to become unstable and drives turbulence to form around the separatrix, which causes a 
stochastisation of the field lines in this area. This process could have significant implications 
for boundary layer theories of tearing stability, as even islands evolving in the presence of 
small gradients can cause the generation of turbulence and stochastisation. 
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